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Abstract 

Quantum random walks have been much studied recently, largely due to their highly nonclassical 
behavior. In this paper, we study one possible route to classical behavior for the discrete quantum 
random walk on the line: the use of multiple quantum "coins" in order to diminish the effects of 
interference between paths. We find solutions to this system in terms of the single coin random 
walk, and compare the asymptotic limit of these solutions to numerical simulations. We find exact 
analytical expressions for the time-dependence of the first two moments, and show that in the long 
time limit the "quantum mechanical" behavior of the one-coin walk persists. We further show that 
this is generic for a very broad class of possible walks, and that this behavior disappears only in 
the limit of a new coin for every step of the walk. 
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I. INTRODUCTION 



In the classical discrete random walk, a particle is located at one of a set of positions. 
In response to a random event (such as the flipping of a coin), the particle moves in one 
direction or another. This sequence is repeated, and the motion of the particle is analyzed 
statistically. These systems provide good models for diffusion and other stochastic processes. 

Considerable work has been done recently on quantum random walks. Quantum walks 
differ from classical walks in that their evolution is unitary and therefore reversible. Two 
classes of models have been studied: continuous [TJ, |3| and discrete ||, [5], ^ ^, [II], [TT|, 
0, 0, 0, 0, [Rj, 0, 0, [HJ [ZD|, |2T], f22j, random walks. This paper is concerned solely with 
the discrete case. 

A naive generalization of the classical random walk behaves in a rather boring way: if 
its evolution is unitary, it can only move in one direction and keep on going that way ||. 
To produce nontrivial behavior, we introduce an extra "coin" degree of freedom (usually 
a single quantum bit) into the system. Just as in the classical random walk, the outcome 
of a "coin flip" determines which way the particle moves; but in the quantum case, both 
the "flip" of the coin and the conditional motion of the particle are represented by unitary 
transformations; thus, there can be interference between different classical paths. 

This paper is concerned with quantum walks on the infinite line. The particle is initially 
located at position x = and is free to travel off to infinity in either direction. We will 
primarily consider the position of the particle after some number of flips t. At this point, 
we assume that the position is measured. We will look at both the probability distribution 
p(x, t) = | (x\ip(t)) | 2 , and at the long-time behavior of the moments (x) and (x 2 ) as functions 
of t. 

Unitary random walks behave quite differently from classical random walks. For a classi- 
cal walk, p(x, t) has approximately the form of a Gaussian distribution (actually, a binomial 
distribution), with a width which spreads like ^/t; the variance x 2 — x 2 grows linearly with 
time. The variance in the quantum walk, by contrast, grows quadratically with time; and 
the distribution p(x, t) has a complicated, oscillatory form. Both of these are effects of 
interference between the possible paths of the particle. 

It should be possible to recover the classical behavior as some kind of limit of the quantum 
system. There are two obvious ways to regain the classical result. If the quantum coin is 
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measured at every step, then the record of the measurement outcomes singles out a particular 
classical path. By averaging over all possible measurement records, one recovers the usual 
classical behavior flTT |. 



Alternatively, rather than re-using the same coin every time, one could replace it with a 
new quantum coin for each flip. After a time t one would have accumulated t coins, all of 
them entangled with the position of the particle. By measuring them, one could reconstruct 
an unique classical path; averaging over the outcomes would once again produce the classical 
result. It is actually unnecessary, however, to measure the coins. Simply by tracing them 
out, one leaves the particle in the mixed state 



P 



^2p(x,t)\x)(x\ , (1.1) 



where p(x, t) is the probability distribution obtained after t steps of a classical random walk. 

It is therefore certainly possible to recover the classical limit. But these two approaches 
give two different routes from quantum to classical. We might increase the number of coins 
used to generate the walk, cycling among M different coins, in the limit using a new coin at 
each step. Or we might weakly measure the coin after each step, reaching the classical limit 
with strong, projective measurements. This is equivalent to having the coin decohere with 
time. 

In this paper we look at quantum random walks with multiple coins. (The case with a 
decoherent coin is considered elsewhere [23, 24 |; systems with decoherence of the particle 



have also been considered [|20| , |22| , P5|.) We begin by reviewing the approaches to solving the 
single-coin walk. We extend these to walks with multiple coins, and derive expressions for 
the amplitudes of an M-coin walk which we can evaluate numerically. We look at the long- 
time behavior of the moments of position both numerically and analytically, and compare 
it to the usual classical random walk. 



II. THE ONE-COIN QUANTUM WALK 

Since our results will utilize the amplitudes for the single- coin quantum walk, we first go 
quickly over the single-coin derivation. We then make use of this in the following section to 
write down exact expressions for the multi-coin case. 

There are two main approaches to analyzing quantum walks. One is to use Fourier 
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analysis, which gives a good qualitative insight into the behavior of the system, but gives 
solutions for the wavefunction in terms of some rather unpleasant integrals. Indeed, no 
exact solution to them is known, and so they must be approximated in the limit as t — > oo. 
The approximation methods become progressively more impractical as the number of coins 
increases; we examine these asymptotic expressions in Appendix B. 

The other approach employs combinatorics. This gives expressions which are opaque, 
but are exact for all times. However, evaluating these expressions exactly is prohibitive for 
long times, and so approximation methods are again required for the long time limit. 



A. Fourier Analysis of Quantum Walks 

Following the analysis in || we consider the Hadamard walk, in which each coin performs 
the evolution 

|R> - £|R) = (|R) + |L» 

|L> „ H\L) = -L (|R) _ |D) (II.l) 

at each time step for which that coin is active, where R and L can be respectively thought of 
as the "heads" and "tails" states of the coin, or equivalently as an internal chirality state of 
the particle. The value of the coin controls the direction in which the particle moves. When 
the coins shows "R" the particle moves right; when it shows "L" the particle moves left. 

Let {\x)} be the position states of the particle, where the values x are integers. We define 
a unitary shift operator 

S\x) = \x + l), S~ l \x) = S ] \x) = \x-l) . {11.2) 

If Vr = \R)(R\ and Vl = \L)(L\ are projectors onto the two states of the coin, then one 
step of the quantum walks is given by the unitary operator 

E=(S®Vr + &®Vl)(I®H) . (11.3) 

If the initial state of the position and coin is \^o), then after t steps of the walk the state is 

|*(t)> = #|tf ) • (n.4) 

For the purposes of this paper, we will usually assume that the coin starts in the state \R). 
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Consider the wave function of the position of the particle. For a one coin walk, this 
wavefunction has just two amplitude components, which are labeled by their chiralities: 

9{x,t)= l aR{x ^) . (II.5) 
\a L (x,t) J 

In Dirac notation this is 

Mt)) = ^2\x)<8>(a R {x,t)\R) + a L (x,t)\L)) . (11.6) 

X 

We would like to solve for aLn(x,t) explicitly. We can do this by diagonalizing the 
operator E. Following the analysis in 0, we use the spatial discrete Fourier transformation, 

5>L,R(k,t) for k E [— 71, 7r] 

/ n dk 
\k) ® (a R (k,t)\R) + a L (k,t)\L)) — , 

oo 

—ikx 



aL,R(k,t) = ^ a L,R{ x ^) e ' 

x=— oo 
oo 

\k) = \*y kx > ( n -7) 

x=— oo 

where the usual integral is replaced by a sum, as the variable x is discrete. We note that 
these "momentum" states \k) are eigenstates of the shift operator S and its inverse: 

S\k) = e~ ik \k) , S ] \k) = e lk \k) . (II.8) 

These states are not normalizable; but we can think of them as the limit of the case where 
the number of position states is large but finite. 

In this momentum basis, the effect of the evolution is given by the matrix 

1 ( e~ ik e~ ik \ 
k y/2 I e ik -e tk J ^ 9 ' ) 

Since this is just a 2-by-2 matrix, it is easily diagonalized, and the final state found in 
terms of its eigenvectors and eigenvalues. The method of solution is now straightforward in 
principle, though complicated in detail. One first represents the initial state in the {\k)} 
basis; applies the diagonal form of the matrix t times; and re-expresses the result in the 
original {\x)} and \R), \L) bases. 
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Nayak and Vishwanath || carried out this program and got exact results. For an initial 
state |* ) = \0)\R), 

t * l + (-l) 2+t r dk / cosA; \ , Ult+kx) 

2 J_ 7T Zn \ Vl + cos 2 kj 

aR ^t) = 1 + { ~ 1)X+t — ; e-"^' . (11.10) 

1 ' J 2 y_ w 2tt Vl + cos 2 fc 1 ^ 



In the exponents, we define the frequencies Uk by sincu^ = y 1/2 sin fc, taking uj^ G 
[-7r/2,7r/2]. 

These integrals, while exact, are difficult to evaluate. Nayak and Vishwanath found 
an approximation to these results, valid in the limit of large t, by making an asymptotic 
expansion of the integrals, to get the approximate expressions 

1 + (-l) x+t 



a L (x,t) 



a R (x,t) 



^8nt(l-x 2 /t 2 ) v / l-2x 2 /¥ 

x(l - x/t) (expi(t(f)(x/t) + tt/4) + c.c.) 
1 + {-l) x+t 



^87rt(l-x 2 /t 2 ) v / l-2x 2 /t 2 
x {{x/t + - 2x 2 /t 2 ) expi(t(f>(x/t) + tt/4) + c.c.) . (11.11) 

The results are conveniently represented in terms of a variable a = x/t. This approximation 
is valid only within the range < a < a/1/2; outside this range the amplitude 

becomes negligibly small. The phase <fi{a) is given by 



'l-2a 2 /l-2a 2 
(a) = arcsin \j - - ^ + a arcsin J - - . (11.12) 



B. Combinatorial Analysis 

We can do a completely different analysis of the one-coin walk, using a combinatorial 
rather than a Fourier argument, similar to that done in || and ||. Suppose we start with 
x = and the coin in state \R). We see how the state evolves after the first few flips of the 
coin: 



l*o> = |0>|J2) , 



|¥(1)> = v/l/2(|l)|i?) + |-l)|L)) , 

|*(2)> = VW(m\R) + |0>|L> + |0)|i?) - | - 2>|L» . (11.13) 



At time t there will be 2* terms, each with amplitude ±-y/l/2*. Each of these terms corre- 
sponds to a possible path for a classical random walk. To get the amplitude for a particular 
position x, conditioned on the coin showing R or L, one must sum up the amplitudes for all 
the paths which end at the position x with the appropriate coin face showing. This amounts 
to a nice problem in combinatorics. 

For the coin to reach position x at time t, it must have moved to the left a total of 
Nl = (t — x)/2 times, and to the right a total of Nr = (t + x)/2 times. For a term 
corresponding to a given classical path, it turns out that the phase depends only on the 
values of Nl and Nr, and the number C of clusters of consecutive L flips. The details of 
the derivation are given in Appendix A, so we just quote the answer here: 



a L (x,t) 
a R (x,t) 



c=i v 



'"N" 

E< 

C=l 



\N L -C 



N L -\ 
C-l 



N R 
C-l 

N R 

c 



(11.14) 



where the summation is to Nl for x > and to Nr + 1 for x < (but only for Nl ^ 0; 
for Nl = the amplitude is always 2~*/ 2 ). We can, of course, re-express Nr and Nl as 
functions of x and t. 

If we had started instead with the coin in the state |L), the amplitudes would then become 



b L (x,t) 
b R (x,t) 



1 



,c=i 

■ u N „ 

E< 

p=i 



1 



\N L -C 



C-l 

N L -l 
C-l 



R 



N R \ N r -2C + 2 
C-l) N, 

Nr\ N r - 2C 
C J Nr 



(11.15) 
is done 



where the summation is to Nl for x > and to + 1 for x < 0, and Nl 
separately, as before. 

The solution given by Eq. ( 11.14 ) looks very different from that given by Eq. ( ILldj) ; it 
is rather extraordinary that they should describe exactly the same outcome. Moreover, a 
long-time approximation similar to Eq. ([Ll.ll|) can be derived from Eq. ( |LL.14 ), by using the 
asymptotics of Jacobi polynomials. See || for details. 

If the coin starts in an arbitrary state 



ikrit = a\R)+l3\L) 



(11.16) 
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the amplitudes will be 



g L (x, t) = aa L (x, t) + (3b L (x, t) (H.17) 
g R (x, t) = aa L (x, t) + f3b L (x, t), (11.18) 



respectively. 



III. WALKS WITH MULTIPLE COINS 



For a walk driven by M coins, the wavefunction has 2 M components; for example, a two 
coin walk will be described by the wavefunction 



\*(x,t)) 



( a RR (x,t)\ 
a RL (x,t) 
a LR (x,t) 

\a LL (x,t)J 



("I 1) 



At each step, the particle moves in the direction dictated by the coin that is active at that 
step, with the other coins remaining inert until it is their turn once again. We will initially 
assume that we cycle through the coins in a deterministic, regular way, although we will 
discover later that we can relax that condition in certain specific ways. 
The unitary transformation which results from flipping the mth coin is 



E m = (S <g> V Rm + & <g> V Lm ){I ® H„ 



(III.2) 



where H m is the Hadamard transformation on the mth coin, and VL,Rm is the projector onto 
the mth coin being in state L or R, respectively. If we cycle among the coins, doing a total 
of t flips (t/M flips with each coin), then the state will be 



\y(Mt)) = (E M ---E 1 ) t / M \y ) . 



(III.3) 



We have simulated this system numerically for different numbers M of coins. The prob- 
ability distributions p(x,t) agree with the classical results only up to t = M; beyond that, 
they diverge sharply, with the multicoin distributions exhibiting highly oscillatory behavior 
and rapid spreading (linear with t), similar to the behavior of the single-coin case. (See 
figure 1.) 
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FIG. 1: These figures plot the particle position distribution p(x) vs. x at time t = 500 for multicoin 
quantum random walks with 1-5 coins. In each of these plots, all the coins begin in state \R), which 
is the reason for the asymmetry of the distributions. Note that only even points x are plotted, 
since p(x) = for odd values of t — x. 
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How can we understand these M-coin results? An important thing to notice is that these 
M unitary transformations {E m } all commute: 

[E m ,E n ]=0. (III.4) 

Because of that, the order in which they are listed is irrelevant. We can equally well write 

|tt(Mf)) = (E M f M ■ ■ • (A)* /M |*o) • (IH.5) 

We can either let the walk evolve for t steps using M coins, or we could use one coin for 
t/M steps, and then replace the coin with another one, and use that for another t/M steps, 
and so on, until we've used all M coins. The final probability distribution will be the same. 

Suppose we decide to use our coins cyclically. Then it can be seen that at the end of the 
walk, the M coins act as a quantum memory of the last M moves that the particle made. 
The "which path" information from before then is lost, so we will see interference between 
paths arising from moves dating from before that time. 

Now consider M walks controlled by one coin each, concatenated in the following recipe: 
Do the first one-coin walk. Measure the final state of the coin only, doing nothing that would 
give any further information about the position of the particle. Reset the coin state to the 
appropriate starting state for the next walk, and use the conditioned state produced by the 
measurement as the starting state for the next one-coin walk. Repeat this procedure until 
you've concatenated M walks in this way. The final probability distribution obtained after 
the last coin has been measured will be exactly the same as the on obtained after measuring 
the M coins from the single cyclic M coin walk. 

How can we use this to solve the multicoin walk? Assume (for the moment) that the 
coins all start in the state \R), and the particle starts in the state |0). We flip the coins t/M 
times each (with t/M even for convenience). After each coin has been flipped t/M times, 
we measure it, and find it in state L or R. We then define conditional evolution operators 
for the particle position: 

t/2M-\ 

A L (t) = a L (2j,t/M)^ 

j=-t/2M 
t/2M 

A R (t) = a R (2j,t/M)S 2j , (III.6) 

j=-t/2M+l 
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where x = 2j. If we flipped M coins t/M times each and then measured them, getting N 
results R and M — N results L, then the particle must have ended up in the (unnormalized) 

state 



N 



A R (t) A L (t) 



M-N 



(III.7) 



The probability of the particle arriving at the final position x at time Mt is then p(x, Mt) 
\^f(x, Mt)\ 2 , which is 

M 



p(x, Mt) = 



N=0 



M 
N 



N 



\(x\ A R (t) A L (t) 



i M-N 



|o)l 2 . 



(III.8) 



Note that this distribution doesn't depend on the order in which the coins were flipped, nor 
the order in which the measurement results occurred, provided that the measurement on 
coin m was made after it had been flipped exactly t/M times. 

For the amplitudes a,L,R(x,t) in (|III.6| ), we can use any of the expressions flll.lOD , (|[I.11D , 



or ( 11.14 ), whichever is convenient. We can also calculate them numerically; Eq. ( III. 8 



provides a much more practical method of numerical calculation for multicoin systems than 
direct simulation, which requires a Hilbert space that grows exponentially in dimension with 
the number of coins. By contrast, the difficulty of evaluating ( |1H.8|) grows only linearly with 
the number of coins. We have used this to simulate systems with large numbers of coins, 
which would otherwise be impractical; see, e.g., figure 2. 

One can find asymptotic expressions for p(x, t) at long times which are analogous to 
the single-coin expressions ( |11.11| ), though the difficulties mount rapidly with the number 
of coins. We plot the approximate solution for two coins in figure 3, which matches the 
numerical results of figure lb quite well. For details of the calculation, see Appendix B. 

While having points of similarity, these multicoin walks differ from the tensor product 



walks studied in ||11|| , in which each component in the tensor product corresponds to a 
different spatial dimension. There are also points of similarity to quantum Parrondo games, 
as studied in 01. 



IV. MOMENTS OF THE DISTRIBUTIONS 



While the probability distributions p(x, t) for the quantum random walk clearly differ 
markedly from their classical counterparts, it is difficult to use this to give a quantitative 
criterion for quantum vs. classical behavior. In particular, we would like to know if the 



11 



0.03 



0.025 



0.02 



0.015 



0.01 



0.005 




200 



FIG. 2: This figure plots the numerical results for p(x,t) at t = 200 for 20 coins, each initially in 
the state \R), using the formula ( |III.8D . As we see, the highly oscillatory form of the probability 
distribution persists even for large numbers of coins. (Only even points are plotted.) 



quantum walks become "more classical" in any meaningful sense as we let M — *■ t. It would 
be useful to look at one or two numbers instead of the entire probability distribution in 
answering this question. 

One obvious criterion suggested by the single-coin case is to look at the long-time be- 
havior of the moments of p(x,t), and see if (for instance) the variance increases linearly or 
quadratically with time. Since it becomes more and more difficult to calculate p(x, t) as we 
increase M, this might seem like an unhelpful approach; but as it happens, we can extract 
expressions for the long time behavior of the moments without having to evaluate p(x, t) 
itself [El. 



Note that in the following two subsections, we fix the number of coins M (or more gen- 
erally the dimension D of the "coin") and then go to the long-time limit. What constitutes 
a long time will depend on the choice of M or D. In the final subsection, we let the number 
of coins used (or the dimension of the "coin" ) increase with time, and look at the long-time 
limit in that case. 
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FIG. 3: This figure plots p(x, t) at t = 500 for the two-coin quantum walk, using the long-time 
approximation derived in Appendix B. Comparison to figure lb shows that this result matches 
direct numerical simulation very closely. 

A. The general case 

Let's consider a very general linear random walk, which includes the multicoin 
a subcase. Let the "coin" be a D-dimensional system with an initial state |$o); let ^Pr,l 
be two orthogonal projectors on the Hilbert space of the coin, such that Vr + Vl = I, and 
TrVn = TtVl — D/2. We also define a unitary transformation U which "flips" the coin. 
Then one step of the quantum random walk is given by the unitary operator 







-400 



-200 








400 



x 




(IV.l) 



where S, are the usual shift operators ( |Tj ) on the particle position. The full initial state 



of the system (particle and coin) is 



(IV.2) 



We can identify the eigenvectors \k) of S, S^, with eigenvalues exp(=pfc), as denned in 
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( p.7|) . In particular, 

io> = r %) . (iv.3) 



In the k basis, the evolution operator becomes 



E{\k)®\®)) = \k) ® (e- ik V R + e ik V L ) U\$) , 

= \k)®U k \$) , (IV.4) 



where U k is a unitary operator on the coin degree of freedom. 

Let the quantum random walk proceed for t steps. Then the state evolves to 

|*o> = |0>® \%) -> ^|^o) 

^ dk,,, /*\* 



2tt 

The probability to reach a point x at time £ is 



\k)®(u k ) \%) = \^(t)). (IV.5) 



p(x,t) = (tf(t)| f J) 

1 /■ „ r , / ,~, ■:■ \ ' / , . \' 



dk / dk'(k\x)(x\k'){® \(ul) fa) |$o) 



(2^) 2 

= ji^I dk I dk ' e ~ ix{k ' k,) (^\ (uiy (u»y i*o> • (iv.6) 

This will, in general, be difficult to calculate. However, all we are interested in are the 
moments of this distribution: 

= (^E^/ dfc / ^'e-^- fe ')($ | (tfj)' (t^)*|$o> • (IV.7) 

We can then invert the order of operations and do the x sum first. This sum can be exactly 
carried out in terms of derivatives of the delta function: 

x m e- ix{k ~ k,) = (-i) m 5 (m \k - k') . (IV.8) 

X 

Inserting this result back into our expression for (x m ) t yields 

(-«? 



(x m ) t = dk dk'5^\k - k'){%\ (ffy (U k ) |$ ) . (IV.9) 
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We can then integrate this by parts to obtain expressions such as 



(x m )t 



-i 



2tt 
2tt 
2tt 



dk j dk'5^\k - k')(%\ (uty (u k ,y i$ ) 



dk / jfe')<$ | (C/J 



dk r 



dk' m 
l$o) 



l*o> 



For the first moment we get 



dk \ k . 



l$o> 



where 



dU k 
dk 

dUl 
dk 
Z 



ie~ lk V R + ie lk V L ) U = -i{V R - V L )U k 



-iZUk 



iui{v R -v L ) = iuiz 



V R -V L = l-2V L . 
Substituting this into ((EV.llQ gives us 



3=1 



We can carry out a similar integration by parts to get the second moment: 



2n 
1 

2n 



l$o> 



i EE / (^) J ^ ^ l$o) • 
i=i i'=i ^ 

Can we now evaluate these expressions? Let us suppose that we can find the eig 
\4>ki) and corresponding eigenvalues exp(i9 k i) of U k . We expand the initial state 



After t steps 



|$o) = E Cm I^ • 



Substituting this into the equation for the first moment ( |IV.13| ) we get 



( x )t 



- j dk^2c* kl c k i f ((f) k i\V L \(f) k i f )^2e 
71 i,v i=i 



i(0kl'- e kl)j 
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If the unitary matrix is nondegenerate, then most of the terms in ( |iV.17| ) will be oscillatory; 



hence, over time, they will average to zero. Only the diagonal terms in the above sum are 
nonoscillatory. We can therefore write 

(x) t = C\t + oscillatory terms , 

Ci = -1 + - / dkJ2\cki\ 2 (<f>ki\V L \<f>ki) ■ (IV.18) 
71 J i 



Making the same substitutions in the equation for the second moment ( |1V.14| ), we get 



(x 2 ) t = y [dkJ2 4^<^%^<^H^I<M E e^'-^e^'- 6 ^' . (IV.19) 

71 1,1', l" j,j'=l 

Once again, most of these terms are oscillatory. There are two sets of nonoscillatory terms: 
terms with I = I' = I" (which give a quadratic dependence on t) and terms with j = f and 
I = /' (which give a linear dependence). Therefore we can write the second moment as 

(x 2 ) t = C^t 2 + oscillatory terms + 0(t) , 

C 2 = 1 - - / dk ^ \Ckl\ 2 {(t)kl\VL\4>kl){4>kl\pR\4>kl) ■ (IV.20) 

Of course, if the spectrum of is degenerate, we will have to modify ( PV.20|) and (|IV.18|) to 
include appropriate cross terms. This does not, however, alter the qualitative behavior. We 
see that generically in the long time limit, the first moment of the quantum random walk 
on the line will undergo a linear drift, and the variance will grow quadratically with time, 
so long as the coin is a finite-dimensional system. 

B. The multicoin model 

Now let us specialize to the case of our multicoin model. In this case our coin is a tensor- 
product of M 2-level coins, with a Hilbert space of 2 dimensions. The "flip operator" is 



U 



®M-1 



(IV.21) 



where H is the usual Hadamard operator ( [II. 1|) and P is a cyclic permutation of the M 
coins: 

P{\^o) ® |Vi> ® • • • ® IV>m-i)) = ® • • • ® IV>m-i) ® |^o)) • (IV.22) 
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The two projectors onto the flip results are 

V L = \L){L\®I mi - 1 . (IV.23) 
When we switch to the k representation, these give us the effective unitary evolution 



P , (IV.24) 



where Hk is given by ( |II.9|) . In the case M = 1, Uk = Hk- 



The eigenvectors of are 

|±) = -f= (l + cos 2 A;± cos kyl + cos 2 k ) I , (IV. 25) 



— e 



with eigenvalues e l( * n+UJk \ e luJk , respectively, where we define cuk to satisfy 

sincufc = —=sink , — ir/2 < uj k < tt/2 . (IV. 26) 

V2 

A basis for the full M-coin space is given by product vectors of |+) and |— ). (E.g., for 

M = 2 the basis vectors would be | ),| — h),| H — ),| + +)•) Can we find the eigenvectors 

and eigenvalues of Uk in terms of this basis? 

Note first that if we apply Uk to such a basis vector, we get another basis vector multiplied 
by a phase. Second, if we apply Uk to a basis vector M times, we get the same basis vector 
back, multiplied by a phase. 

Two obvious eigenvectors appear at once: the vectors 

I + + = |+)^ M , 

| ) = \-f M , (IV.27) 

are both eigenvectors of Uk with eigenvalues e^ n+UJk \ e~ tuJk , respectively. Let us suppose for 
the moment that M is a prime number. Then for all basis vectors |0) other than those two, 
the vectors 

Uk) 3 \4>), 0<j<M-l (IV.28) 

are all distinct. Together they span an M-dimensional subspace which is preserved under 
the action of Uk- If \4>) contains exactly m |— )s and (M — m) |+)s, then every vector in this 
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subspace is an eigenvector of (Uk) M with eigenvalue expi((M — m)(iT + Uk) — rau^). This 
implies that the subspace is spanned by M eigenvectors of Uk with eigenvalues 

A n = e i(-^ fc +(M-m)( 7 r+a; fc )+2 7 rn)/M ^ < n < M . (IV.29) 

Let's restrict ourselves to this M- dimensional subspace for the moment. Starting with 
our original basis state \<p), we label the basis states which span this space 

|&) = P>|0). (IV.30) 
The eigenvectors of U^ must have the form 

M-l 
V 3=0 

Plugging the above expression and equation ([[V.29|) into the equation 



) (IV.32) 

gives the result 

Vni+i — v nj - {muj k — (M — m)(7T + uJk) - 27m) /M + 8 j+ i , (IV.33) 

where we define 9j to be e~ lulk if the first factor of \4>j) is |— ) and to be e^ 7r+aJfe ^ ) if the first 
factor is |+). Since the overall phase is arbitrary, we take u = 0, which gives us 

3 

Vnj = -jj ( muJ k - (M - m) (tt + u k ) + 27m) + d f ■ {IV M) 

f=i 

For each value of m, 1 < m < M — 1, there are ) such M-dimensional subspaces 
which are preserved by the action of Uk', each eigenvalue is thus -g(^)-fold degenerate. 
However, this doesn't complicate the expressions (|IV.i8| ) and ( |IV.20| ), because the cross- 
terms vanish: 

{<t>ki\v R \M = (failAlM = o (iv.35) 

for all such degenerate eigenvectors \(f>ki), \4>kv)- Thus, we can use the expressions we've 
already derived. 
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Using the results in Eq. ( |IV.34|) we can calculate the necessary matrix elements 

1 



(Xn\V L \Xn) 



M 
1 

M 



m (-\L){L\-) + (M-m){+\L){L\+)} 
ra(\/l + cos 2 k — cos k) 2 



2(1 + cos 2 k — cos kyl + cos 2 k) 

(M - m)( v / l + cos 2 A; + cos kf 
2(1 + cos 2 k + cos k\Jl + cos 2 k) 



{Xn\P R \Xn) 



M 

1 

M 



[m(-\R)(R\-) + (M — m)(+\R)(R\+)} 



m 



2(1 + cos 2 k — cos fcvl + cos 2 k) 
M -m 



2(1 + cos 2 k + cos kyl + cos 2 fc) 
We also need the amplitudes of the coin's initial state. If 



then 



Note that 



|$0> = KM ® • • " ® |-0O> = l^o) 



l($o|Xn)| 2 = |(^o|-)| 2m |(^o|+)| 2(M - m) E -e^-^') 



M-l 



^ e i("»i-"ni') = M^v . 



n=0 



We now have everything we need to get the moments! 



To evaluate the expression (|IV.18|) to get C\ we calculate 



E 



M_1 1 fM\ m ~ L 

VM = M\m) S \(®0\<Pn)\ 2 (<j)n\V L \<Pn) 

m=l ^ ' n=0 

H-K^ol— >| 2 < — |X,> <X,|— > + |^o|+)| 2 (+|^)(^|+) 

E(^)i(^oi-)i 2m i^oi+)r (M - m) 



M-l 



Ckl {0kl 



m=0 



X- 



M 



m{\Jl + cos 2 — cos fc) 5 



2(1 + cos 2 /c — cos + cos 2 k) 

(M - m)( v / l + cos 2 A; + cos/c) 2 
2(1 + cos 2 A; + cos kyl + cos 2 fe) 
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We can simplify this considerably by noting that this equation ( |IV.40| ) has the form 



J2 ( )P m ^ -p) M ~ m [(™/M)(A- B) + B] = pA + (1 - p)B , (IV.41) 



where 



(y/1 + cos 2 k — cos k) 2 



A 
B 



2(1 + cos 2 k — cos ky/1 + cos 2 k) 

(y/1 + cos 2 k + cos k) 2 
2(1 + cos 2 k + cos kyl + cos 2 k) 

and we've invoked some combinatorial identities 

M 



(IV.42) 



m=0 v 7 



M-m 



(P+(1- P )) M = 1 



M 



M 



m 



p m (l -p) M ^ m (m 



m=0 x 7 m=0 

This leaves us with the simple expression 



M_1 //If — 1 \ 

/M) = p ( P m (l -p)^" 1 - 1 = p . {IV A3) 



l(^o|-)| 2 (v / T+~cos 2 £; - coskf 



2(1 + cos 2 k — cos fcyl + cos 2 k) 
K^o|+)| 2 (v / l + cos 2 A; + cosfc) 2 



2(1 + cos 2 k + cos ky/1 + cos 2 fc) 

in which all M dependence is gone! 

Let us choose the initial condition |-0o) = \R)- Then 

1 



(IV.44) 



2(1 + cos 2 k — cos ky/l + cos 2 k) 



(IV.45) 



2(1 + cos 2 k + cos ky/1 + cos 2 k) 

We can plug this into the above equation ( |IV.44|) ; after a bit of algebra, this boils down to 
the result 



2J \Ckl\ 2 (^kl\VL\(pkl) 



2(1 + cos 2 k) 



(IV.46) 



This gives us the final equation for C\ 

1 



Ci 



1 



27T / 1 + COS 2 



-1 + 1/V2, 



(IV.47) 
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FIG. 4: This plots the first moment of position vs. time for multicoin quantum random walks 
with M =1-5 coins, using both direct numerical simulation and the analytical result flTV.47 ). As 
calculated, the long time behavior is essentially independent of the number of coins, and grows 
linearly with time. All coins begin in the initial state \R). 

which exactly matches the observed numerical results. (See figure 4.) 

We can do a similar derivation for the second moment, in order to get C^. For this we 
use 



I Cfcj | 2 (4>m I Vl I <Pki) (4>ki I Vr I <Pki) 

I 

= E( M )i^oi-)ri^oi+)i 2(M - ) 

m=0 ^ ' 

m{yl + cos 2 k — cos k) 2 



M 



2(1 + cos 2 k — cos k\/l + cos 2 k) 

(M - m){y/l + cos 2 k + cos/c) 2 
2(1 + cos 2 k + cos ky/1 + cos 2 k) 



M 



m 



2(1 + cos 2 k — cos fcyl + cos 2 k) 

(M - m) 
2(1 + cos 2 k + cos k\/l + cos 2 k) 



(IV.48) 
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We can again simplify by observing that ( |LV.48|) has the form 

M 



^ (m) 



p m {l - p) M - m [(m/M)(A -B) + B][A+ [m/M)(B - A)} 

M — 1 



m=0 

A/T _ 1 



= AB+p(l-p) — (A — B) . (IV.49) 
If we specialize once more to = |0), we can make similar algebraic simplifications to get 



Co = 1 - — / dk 



1 + 2 cos k 1 cos k 



2tt J 1(1 + cos 2 k f M (1 + cos 2 A;) 2 
5 +^tV (IV.50) 



4^ M Ay/2 ' 

Of course, numerically we calculated not the second moment but the variance of x. Given 
our results, 

(x 2 ) t - (x) t 2 = (C 2 - Cl)t 2 + 0(t) + osc. terms , (IV.51) 

we get the result 

Note that this doesn't vanish as M becomes large! 

The derivation used above was only strictly valid for prime M. However, it turns out 
that the additional complications arising from composite M all cancel out, yielding the 
same expressions for C\ and C 2 . Upon comparison to the numerical results, we see that 
these analytical expressions for the first moment and the variance match the simulations 
extremely well. (See figure 5.) 



C. Constant flips per coin 

We see that in the long time limit, an M-coin quantum random walk with fixed M has 
qualitative behavior similar to that one of the one-coin walk, and is markedly nonclassical. 
A reasonable question to ask is, how long is a long time? Suppose we flip each coin at most 
d times, and add more coins as we go to long times, so that M and t are both growing in 
a fixed ratio t/M — > d. Does this still behave nonclassically? The fact that C% — C\ does 
not vanish as M — > oo makes this conjecture plausible. If so, how big must d be for this 
nonclassical behavior to manifest itself? 

Let us now consider a quantum random walk on a line driven by a sequence of two-level 
coins, each flipped d times. By the reordering principle described in section III, this is the 
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FIG. 5: This plots the variance of position vs. time for multicoin quantum random walks with 
M =1-5 coins, using both direct numerical simulation and the analytical result ( |1V.52| ). All coins 
begin in the initial state \R). 



same as having a single coin which is measured and reset to the initial state |$o) after every 
d steps. We describe this evolution by a superoperator £ which acts on the density matrix 
of the particle and coin: 



p^p' = Cp 



(IV.53) 



where E is the usual unitary evolution (|II.3| ), and 



A = I®\%){R\ , 
Ax = I® \%)(L\ . 

The superoperator C represents d steps of the walk. Note that 



Cp = Tr coin (i f ) d | <g> |$q)($o| 



(IV.54) 



(IV.55) 



We rewrite this in terms of the eigenvectors \k) of S, . In the k basis, a general density 
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operator for the joint particle/coin system is written 



dk f dk' 

P= I 7T I tHM*'!®***, (IV.56) 



2tt J 2tt 

and the evolution superoperator becomes 

n=R,L 

= J^J ^\k)(k'\®£ kk 'Xkk> , (IV.57) 

where Ckk' is now a superoperator on the coin degree of freedom alone. 
The initial state of the system (particle and coin) is 

Po = l^oX^ol = |0)(0| <g> |$o)($o| 

= j^j d ^\k){k>\®\%){%\. (IV.58) 

Let the quantum random walk proceed for t steps. Then the state evolves to 

Pt = J^J d ^\k)(k^Ci d ,\%)(%\. (IV.59) 
The probability to reach a point x at time t is 
p(x,t) = Tr{(\x)(x\®I)pt} 

= J^fI dk I ^(k\x)(x\^{c^ )(^ \} 

= (2^/**/ ^' e "^ fc ' )Tr {4 / ,^o)(<fo|} • (IV.60) 
We are interested in the moments of this distribution. 

(x m ) t = £*>0m) 

X 

= J^yY, xm j dk j ^'e-( fe - fc ')Tr{4^|$o)($o|} ■ (IV.61) 
Just as before, we invert the order of operations and do the x sum first, which yields 

(n t = ^ jdkj dk'5^\k-k')Tr{c^\%)(%\} . (IV.62) 
We can then integrate this by parts. In carrying out this integration by parts, we will need 

n \j=0 J 

= -^o)(*o|Tr | {^{H k yZ{Hl)^j (H h )*p(& k )^ , (IV.63) 
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where H k is given by ( |II.9|) and Z by (|IV.12| ). 

Using this, we can carry out the integration by parts for the first moment to get 

t/d-l 



i- ]T ! dkTr {Z kd 4 k \ $ )($o 
1 /" 

— J dkTr{z kd \%)(^ 



where we define 



j=0 

t f dk , . - . 



1=1 



(IV.64) 



(IV.65) 



Generically, then, we see from (|IV.64j) that the first moment of position grows linearly with 
time provided (^ \Z kd \^ ) is nonzero. 

We can carry out a similar integration by parts to get the second moment: 

/'t/d-l (2 1 
dk J2 Tr { (Z*d) 4J*o><*o|)f 

. 3=0 L } 

t/d-l j-l 

+ E E^{ i ^i/((4>o)(*oi)4d+4 d 4' fe i^o>($Qi) 

J'=0 j'=0 

(t/d)($o| (Zkd) 2 |$o) + (t 2 /^ 2 - t/d)(%\Z kd \<!> } 2 . (IV.66) 



— / dk 

2vr 



So here we expect the second moment to grow quadratically with time, qualitatively like the 
single-coin walk. 

Let's choose the initial condition |$ ) = \R), an d look at a few values of d. We are 
interested in the coefficients C\ and C2, where 



(x) t = dt , 
(x 2 ) t = C 2 t 2 + 0(t) . 



(IV.67) 



From the equations ( |1V.64| ), ( |i V.66| ) , we get 

„ 1 f dk 

L\ — 



d / 2tt 



R\Zkd\R) 



c 2 = U "wuv ■ 



(IV.68) 
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From the definitions of Z^d and we can work out the values of C\p. in specific cases. 

For d = 1, the matrix element (R\Zki\R) = 0, so both C\ and C2 vanish. We recover the 
classical case. 

For d = 2, the matrix element (R\Zk2\R) = cos2A;, so C\ vanishes in this well. 
However, C2 = 1/8 7^ 0. So even with only two flips per coin, the variance already grows 
quadratically with time. 

For d = 3, the matrix element [R\Z^\R) = cos2A; + sin 2 2 A;, so both C\ and C2 are 
nonzero. In this case, C\ = — 1/6 and C2 = 7/72. 

If we let both d and t/d become large, then we can evaluate ([1V.64J) and ( |1V.66| ) to get 
C x -»• -1 + and C 2 -»• 1 - 5/^32, which are in agreement with our earlier results for 

large M in the multicoin case. 

V. CONCLUSIONS 

We have examined one possible path from quantum to classical behavior: the use of 
multiple coins (or, more generally, higher-dimensional systems) to drive the walk and reduce 
the effects of interference. We have seen that quantum behavior, as typified by quadratic 
growth of the variance with time, persists except in the extreme limit of a new coin for every 
step. Furthermore, quadratic growth of the variance seems to be a generic feature of such 
unitary walks. 

From the multicoin example we might speculate that classical behavior is only recovered 
in the limit where the coin system retains enough information about the walk to reconstruct 
a unique classical path for the particle. Since there are 2* such paths up to time t, the coin 
must have a Hilbert space dimension which grows exponentially in time in order to exhibit 
classical behavior. 

Another plausible route to classical behavior adds decoherence to the coin, which also 
effectively suppresses interference effects. The results of this study are presented elsewhere 

& 
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APPENDIX A: COMBINATORIAL DERIVATION OF ONE-COIN WALK 

Assume that the coin starts at x — 0. After t coin flips, the particle might have reached 
locations anywhere between x = —t, — t + 2, . . . , t — 2, t. The number of paths that finish at 
a point x after t steps is 

#PathS =(( t -l)/2)' (A1) 
Paths ending with the coin in state \L) cannot interfere with paths ending with the coin in 

state \R). 

A given path can be written as a sequence S of t symbols, L or R like this one: 

L RRL RL_L_LRR... (A.2) 

a cluster 
of 3 Ls 

We let N L (S) be the number of L's in S and N R (S) the number of i?'s. Then 

N L (S) + N R (S)=t, N R (S) - N L (S) = x , (A.3) 

so we can write 

N R (S) =(t + x)/2, N L (S) = (t- x)/2 . (AA) 

Each path has an associated phase factor, which can take the values ±1. A factor of — 1 
is acquired whenever the coin flips two successive lefts, so a pair LL contributes a phase 
factor of —1, the sequence LLL produces a factor of +1, and LLLL results in a —1, and so 
on. Each cluster of m Ls contributes a factor of (— l)™" 1 to the phase for that path. 

Let C(S) be the number of L-clusters in S. (An L-cluster must contain at least one L.) 
Swapping an L from one cluster to another changes the phase contributed by each cluster 
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by —1, thus having no overall effect on the phase of the path as a whole. Therefore, only 
the number of L-clusters and the total number of Ls matters. 

If the number of L-clusters is C and the total number of Ls is N L , then these could be 
arranged thus 

L L ... L LL . .. L . (A.5) 

C-l singletons N L -{C-1) 
consecutively 

This path will have a phase of 

<p(S) = (-i)^-(c-D-i = ( A . 6 ) 

So when Nl — C is odd the overall phase is —1. 

Each cluster must have at least one R between itself and the next cluster. Call these 
groups of Rs "partitions." Every sequence S is then a succession of alternating clusters and 
partitions. Let 

P(S) = # of partitions in S. (A.7) 

Clearly C(S) cannot be higher than Nl(S), and P(S) cannot be higher than Nr(S). Also, 
it must be the case that 

P(S) = C{S) or C{S) ± 1. (A.8) 

Allocating Nl Ls among C clusters can be done in ( i ^f_T 1 1 ) ways; similarly, Nr Rs can be 
allocated among P partitions in ( N pZi) ways. For a given C and P, then, the total number 
of paths must be 

"-■WW' 

Let us now fix t and x, and hence also N L and N R . What values of C and P are possible, 
and how many paths have each value? 

Case I: x > =^ N R > N L . The number of clusters, C, can range from (only when 
x — t) to N L . The number of partitions P can be C, C + 1 or C — 1, except when C — 0, 1. 

Case II: x < =^ iV L > A"/?. This time C can range from 1 to N R + 1. Likewise, P 
can take values of C, C + 1, C — 1 except at the extremes. 

Case III: x = Nl = N R and C and P can be anything from 1 to Nl = N R . 

Consider now the arrangements of clusters and partitions. There are four distinct possible 
arrangements. 
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Arrangement 1: P = C — 1. 



[cluster] partition [cluster] partition . . . [cluster] 



Arrangement 2: P = C, starting with a partition. This looks like 



partition [cluster] partition . . . [cluster] 



Arrangement 3: P = C, starting with a cluster. 



[cluster] partition [cluster] . . . partition. 



Arrangement 4: P = C + 1. 



partition [cluster] partition [cluster] . . . partition. 



(A.10) 



(A.ll) 



(A.12) 



(A.13) 



We adopt the convention 



0, if b > a, b < 0, or a < 0. 



(A. 14) 



Then we can write that the amplitude to reach x after t coin flips, ending in an L is: 



a L (x,t) 



C=l v 



N R - 1\ (N R - 1 
C-2 J [ C-l 



(A.15) 



where the summation is to Nl for x > and to Nr + 1 for x < 0. 
Likewise, the amplitude to reach x in t flips, ending in an R is 



a R (x,t) 



"N" 



,c=i 



C-l 



iV fi - 1\ /iV fl - 1 

C-l J v c 



(A.16) 



but only if Nl ^ 0; that case (which corresponds to t = x) always has amplitude 2 */ 2 . 
We can use the binomial relations 



N R - 1\ /iVfl - 1 

c-2y \c-i 

- 1\ - 1 

c-i 7 I c 



C-l 

c 



(A.17) 



to simplify these somewhat to give the final equations (p.!4|) . 

If we had started instead with the coin in the state \L), all sequences which begin with 
an L (i.e., arrangements 1 and 3) would pick up an extra minus sign. The amplitudes then 
become those given by the equations ( |II.15| ). 
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APPENDIX B: ASYMPTOTIC APPROXIMATION TO THE FOURIER INTE- 
GRALS 



Following the analysis in , the Fourier analysis of the many coin quantum walk can be 
completed in the same way. We represent the state of the particle in the {\k)} basis, and 
make use of the eigenvectors and eigenvalues of Hk, as given by ( JIV .25 ) and ( |lV.26p . We 
will also use the identity 

±V2e Ti " k - e~ lk = ± v /l + cos 2 {k) - cos{k) (B.l) 

We can now use the tensor-product structure of the evolution matrix and the integral sim- 
plifications in || to simplify the integrals for clr and a^. The expression for a# becomes 

_ cos(A;) \ el{{uk+ ^ kx)dk = +t r ( co S (fc') \ e _ i{Uk , t+klx)dk , 

y/l + cos 2 (A;) J J- n \ v/l + cos 2 (A;') J 

(B.2) 

where k' = k — 7T, Uk> = —Uk and we are free to move the limits of integration because the 
integral is over a whole period of the function. The two terms in the integral for cll can be 
simplified in a similar way: 



_ e iki 



'{y/1 + COS 2 (fc) + COs(fc)) pi{wk+n ) t -kx dk 
((1 + cos 2 (fc)) + cos(fc) y/1 + COS 2 (fc)) 

6 e -i{u k ,t+k' X ) dk , _ /g 3 N 

2y/l + cos 2 (A;') 

Before we can write down a closed form for the M-coin momentum wave-function, we 
need one more construction. This is a timelike parameter that we'll call |r|. Let tj be 
the total number of turns that coin i has. The tensor product structure of the evolution 
matrix for the multicoin system means that the eigenvalues for the combined system are 
just products of the eigenvalues for the one coin system. So we'll see things like 

^ _ e ~iuJ k t-i e i(uj k +Tr)t2 e -iuJ k t Se i(u) k +Tr)t4 e i{u>k+Tr)t5 ^ 



If we now write 



(h -t 2 + t 3 -t A -t 5 ...)= "r" (B.6) 
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we can see that, up to a minus sign, we will obtain expressions of much the same form as in 
the one coin case, but parametrized by r instead of t, where r is defined 



re{±t 1 ±t 2 ±...±t M }, (B.7) 

with each ± independent. The modulus sign arises naturally when the integral expressions 
for the wavefunction are simplified in a way analogous to the one coin case, as r can be 
chosen to be positive without loss of generality. When performing the simplification, it is 
also convenient to introduce the integer j T , which is useful because for each r we get terms 
like 

9/ \ #(tiS ending in L) / \ #(—tiS ending in L) 

- 1 ; cos i h ; os {k) ) (b.8) 

yjl + COS 2 (A;) J \ y/l+COS 2 (k) J 

These are of the form (a + b) r \a — b) s , and j T counts the (a 2 — fe 2 )l r_s l terms in the expan- 
sion of that product. These are of the form ( 1+co 1 s a ( - fc - ) )^~' s ^ and so they only appear in the 
denominator. The term raised to the power of n — 2j T is just the two sets of left-overs from 
that (one multiplied by e - l ( aJ fcM+ fc:c ) anc i the other by ^((^k+^M-kx)^ w j 1 i c } 1 j^yg then been 
combined using the integral simplification identities as before. 

We can now write down the wavefunction for a particle that started out with all the coins 
in the state \R) and ended in them in (any) one of the ( M ) components with n coins in the 
state R : 



\lpn,(M-n)(x,t)) 

E 



2 M 2tt 



1 + (_l)«+|r| rn dk ( e ^(A/-n) ( cos (fc) + ^1 + COS 2 (fc)^ " ^ 

(l + cos 2 (A;)) M / 2 



e -i( Wfc |r[+te) ^ ^ B _g) 



where the sum over r indicates a sum over all 2 M possible signs in ( p.7| ). 

To write down the wavefunction for a particle that started with its coins in any state 
in the computational basis for the coins, we need yet another definition. We will introduce 
the notation o to be the L-sign weight of r, which will be the number of minuses in that 
particular r. So for 

T = tl + t 2 - t 3 + U - t 5 - t 6 + t 7 , 

initial coin state = RLRLLLR, (B.10) 
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we define a = 2. In other words, for all the coins that started out in the state L, count the 
number of minuses in front of the number of turns taken by those coins for that r. Ignore 
any minuses for coins that start in the state R. Without loss of generality, a can be taken 
to be in [0, [-W/2J], as only |r| actually matters in the wavefunction. 

Then the wavefunction for a walk starting with its M coins in a state in the computational 
basis of the coins, where there are initially q coins in the state L and M — q in the state R 
will be 



|V ; n,(M-n):coins(3 ; ) t)) — ^ ] 

T 

" n dk ,„ 



x 



2tt 



1 + 

1 + cos 2 (/c) — cos(fc) 



M-q-cr 



1 + cos 2 (/c) + cos(fc) 



X 



- ikn 



cos(fc) + a/1 + COS 2 (fc) 



n—2j T 



,jikq 



\ 



+ COS 2 (fc)) M / 2 



-i(cOk |t|+A;x) 



fB.ll^ 



/ 



The wavefunction for a general starting state can be obtained from those for the compu- 
tational basis states by forming the corresponding linear combination. 

If r is an asymptotic parameter: i.e., r — > oo as t — > oo then these integrals may be 
approximated using the method of stationary phase. If r = or r = constant (which can 
happen if the coins are used cyclically: for example if we have two coins, then we will get 
terms where t\ — t 2 = or 1 for all time) then the integrals must be performed by some 
other method. It is these constant-r cases that produce the "stationary" central spikes seen 
in the simulations for an even number of coins (more on these below). 

If we restrict to the special case where ti = t 2 = • • • = t m = t/M is an integer, then 
equation (|B.9|) simplifies somewhat to 



dk 

2^ 



LM/2J Ln/2J 

E E 

r=0 j=0 
( 



M\ 1 + (_l)aH-(Af-2r)t/M 



)M 



X 



COs(fc) + y/1 + COS 2 (AO) e iHM-n) 



n-2j 



-n)\ 



\ 



[1 + cos 2 (fc)) M / 2 



-i{oj k (M-2r)t/M+kx) 



(B.12) 
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All of the time-dependent integrals are of the form 

dk 



2tt 



i(/3(fc,O r )|T| 



(B.13) 
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i.e., these are one-parameter families of generalized Fourier integrals, with parameter a T = 
x/t. They can therefore be approximated in the limit as r — > oo by the method of stationary 
phase, as used in ||, and references therein, principally [f26[] and ^7|. When each coin is 



used the same number of times, it can be seen that the spikes are equally spaced between 
-t/y/2 and t/y/2. It can also be shown that the height of the peaks scales like 1/ 

The asymptotic behavior of the time-dependent terms in the above components is as 
follows. Just as in the one coin case, the integral displays three different types of behavior, 
depending on a T , and we refer readers interested in the details of the analysis to ||. For 
the purposes of calculating the moments of the distribution, we can treat the support of the 
integral as confined to the region a T G ( — ^75) to a good approximation, and so we will 
only do the calculations for a T in this range. All the integrals have the same form, and differ 
only by a linear rescaling. 

The values of k corresponding to the two stationary points of ip(oc T ) are ±k aT , where k 
is a function of a T , and 

— (X 

cos(£w) = ~7==- (B.14) 

We also write 

(p(a T ) = (uj kaT + a T k ar ) u kaT = arcsin ( J ^ ) , (B.15) 

and note that 

1 " 1 

I rU(XT 1 

cos(fc) 

1 + V '= = 1 - a T (B.17) 




(l-t#Vl-2a£ (B.16) 



1 + cos 2 (A;) 



e 



ik 



-Or + iy/l - 2a 2 T . (B.18) 



y/\ + COS 2 (k) 

Using these, we can now write down our asymptotic approximations for the time- 
dependent components. Let us consider the two-coin walk as an example. The time- 
dependent components i[)(x,t) for the two-coin wavefunction are then 

\4>rr(x, £)) = ; 1 + (1 _ a y cos (^( a ) t + ^ (B.19) 

27rt(l - a 2 )Vl - 2a 2 
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1 4. r_-n(«+i)* 

\4>lr) = \^rl) 



X 



(a(a - 1) cos(ip(a)t + 7r/4) + (a - 1) Vl - 2a 2 sm(ip(a)t + 7r/4)) , (B.20) 
and 

1 + 



y^27rt(l -a 2 )Vl -2a 2 

x ((3a 2 - 1) cos(y?(a)t + tt/4) + 2oVl - 2a 2 sin(y>(a)i + tt/4)) , (B.21) 



where we have dropped the r subscript on a because there is only one set of asymptotic 
integrals for this system. These expressions can be used to calculate p(x, t) in the long time 
limit; we plot this result in figure 3. 

More generally, for those components which can be obtained using the method of sta- 
tionary phase, we can write 



, . , „ v^ (l + (-ir + '-')(l-a r )("- 2 »(l 

|YV(M-n)(Z,tl, • • • ,*Mj) = / — 

2 M - 1 y / 27r|r|(l-a 2 ) v /r^2^ 



a 2j(M+j T -n/2) 



x cos(v?(a r )|r| + (M — n)k ar + tt/4), (B.22) 

where \ip n ,(M-n)(x, h, . . . , tu)) denotes any of the (^) components of the wavefunction which 
end up with n coins in the state R and (M — n) coins in the state L. 

If M is an even number, the wavefunction will have components tp(x) that are central 
spikes which perform periodic motion for a cyclic walk as the walk cycles through the coins. 
If the coins are used in some other order, these spikes are the final positions of part of the 
wavefunction. In any event, these central spikes make no contribution to the moments of 
the distribution. 

Let us write the total wavefunction if)(x,t) = ip(x) + ip(x,t), where i/)(x,t) denotes the 
rest of the components. For the two-coin case, the central spike components are as follows: 

f-lVi + (-1Y 2 r e~ ikx 

fcW> = (-D h + (-D- f L^W dk + 

8n J_ w 1 + cos 2 (fc) 



^ L_bL / rfA: (B.24) 

8vr J_„ v /1 + cos 2 (A;) 
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f_1 \ti , (_-l)t2 -ik(x-2) 

= LJ_±L_L £ ____<«. (B.25) 

These integrals shrink to zero very rapidly away from x — 0, but they are still visibly 
non-zero out as far as x = ±6 or so. 

These time-independent stationary spike integrals are similar in form to the time- 
dependent ones, but this time the method of stationary phase isn't helpful. That asymptotic 
approximation method can only tell us about their behavior as x — > oo. Since we already 
know that they're independent of time, we can say that these represent something that stays 
in the vicinity of the origin and doesn't go anywhere: so we already know they'll tend to 
zero very rapidly for large x. This means that we're only interested in their behavior for 
a few points in the vicinity of the origin. Since the kernel of the integrand oscillates fairly 
slowly for small values of x, it is practical to evaluate these terms numerically. 
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